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Abstract: We critically reexamine the standard applications of the method of kinematical 
endpoints for sparticle mass determination. We consider the typical decay chain in super- 
symmetry (SUSY) g — > X2 ^ ^ ^ Xii which yields a jet j, and two leptons £j and £j. The 
conventional approaches use the upper kinematical endpoints of the individual distributions 
rriju, mj^(j_o) = nim{mji^,mji^.} and mjK^fii-^ = max{mji^,mj£^}, all three of which suffer from 
parameter space region ambiguities and may lead to multiple solutions for the SUSY mass 
spectrum. In contrast, we do not use rriju, mjn^ig-^ and rrijuj^i^, and instead propose a new set 
of (infinitely many) variables whose upper kinematic endpoints exhibit reduced sensitivity to 
the parameter space region. We then outline an alternative, much simplified procedure for 
obtaining the SUSY mass spectrum. In particular, we show that the four endpoints observed 
in the three distributions m^^, rn^j^^ U fn^ji^ and rn^j^^ + "^^^^ are sufficient to completely pin 
down the squark mass rriq and the two neutralino masses m^o and m^o , leaving only a discrete 
2-fold ambiguity for the slepton mass m^. This remaining ambiguity can be easily resolved in 
a number of different ways: for example, by a single additional measurement of the kinematic 
endpoint of any one out of the many remaining 1-dimensional distributions at our disposal, 
or by exploring the correlations in the 2-dimensional distribution of fn^jf^^ U fn^jn^ versus m^^. 
We illustrate our method with two examples: the LMl and LM6 CMS study points. An ad- 
ditional advantage of our method is the expected improvement in the accuracy of the SUSY 
mass determination, due to the multitude and variety of available measurements. 



Keywords: Supersymmetry Phenomenology. 



Contents 



1. Introduction 2 

1.1 Near-far lepton ambiguity 3 

1.2 Insufficient number of measurements. 5 

1.3 Parameter space region ambiguity 6 

1.4 Posing the problem 7 

2. New variables 8 

2.1 The union rrij^^ U m'j^^ 8 

2.2 The product mj^^ x niji^ 9 

2.3 The sums m'j"^ + rnj'^^ 9 

2.3.1 Kinematic endpoints of "^j£(s)(ct) with a > 1 9 

2.3.2 Kinematic endpoints of 'm'^£(^g^{o() with a < 1 and a 7^ 10 

2.4 The difference - ^i^^^^^l 11 

3. Theoretical analysis 12 

3.1 Our method and the solution for the mass spectrum 12 

3.2 Disambiguation of the two solutions for uib 13 

3.2.1 Invariant mass endpoint method 13 

3.2.2 Invariant mass correlations 14 

3.2.3 Mt2 endpoint method 16 

4. Numerical examples 17 

4.1 Mass measurements at points LMl and LM6 17 

4.2 Eliminating the fake solution for uib 21 

5. Summary and conclusions 24 

A. Appendix: Analytical expressions for the shapes of the invariant mass 

distributions 27 

A.l Dilepton mass distribution mj^ 27 
A. 2 Combined jet-lepton mass distribution 'yn'^jHu) 

A. 3 Distribution of the sum m|^|-^.^(a = 1) 28 

A. 4 Distribution of the difference m^^j^^j (a = 1) 29 

A. 5 Distribution of the product m%i^\ 31 



- 1 - 



j in if 




Figure 1: The typical cascade decay chain under consideration in this paper. Here D, C, B and 
A are new BSM particles, while the corresponding SM decay products are: a QCD jet j, a "near" 
lepton £^ and a "far" lepton £j . This chain is quite common in SUSY, with the identification D = q, 
C = X21 B ^ i and A = Xij where g is a squark, ^ is a slepton, and Xi (X2) the first (second) 
lightest neutralino. In what follows we shall quote our results in terms of the D mass m^) and the 
three dimensionless squared mass ratios Rcd, Rbc and Rab defined in eq. (1.6). 

1. Introduction 

SUSY is a primary target of the LHC searches for new physics beyond the Standard Model 
(BSM). In SUSY models with conserved iZ-parity the superpartners are produced in pairs 
and each one decays through a cascade decay chain down to the lightest superpartner (LSP). 
If the LSP is the lightest neutralino y^, it escapes detection, making it rather difficult to 
reconstruct directly the preceding superpartners and thus measure their masses and spins. 
In recognition of this fact, in recent years there has been an increased interest in developing 
new techniques for mass [1-49] and spin [50-76] measurements in such SUSY-like missing 
energy events. 

Roughly speaking, there are three basic types of mass determination methods in SUSY^. 
In this paper we concentrate on the classic method of kinematical endpoints [1]. Following 
the previous SUSY studies, for illustration of our results we shall use the generic decay chain 
D ^ jC ^ ji^B ji^ijA shown in Fig. 1. Here D, C, B and A are new BSM particles 
with masses rriD, mc, ms and niA- Their corresponding SM decay products are: a QCD jet 
j, a "near" lepton £^ and a "far" lepton £j. This decay chain is quite common in SUSY, 
with the identification D = q, C = B = i and A = Xi, where is a squark, ^ is a slepton, 
and Xi (X2) ^^^^ (second) lightest neutralino. However, our analysis is not limited to 

SUSY only, since the chain in Fig. 1 also appears in other BSM scenarios, e.g. Universal 
Extra Dimensions [77]. For concreteness, we shall assume that all three decays exhibited in 
Fig. 1 are two-body, i.e. we shall consider the mass hierarchy 

mo > mc > mB > mA > 0. (1.1) 
^For a recent study representative of each method, see Refs. [43,47,49]. 
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This presents the most challenging case, in which one has to determine all four masses niD, 

ruB and m,A- 

The idea of the kinematic endpoint method is very simple. Given the SM decay products 
j, in and £f exhibited in Fig. 1, form the invariant mass^ of every possible combination, me£, 
rriji^, rriji^, and mj^^, plot the resulting distributions and measure the corresponding upper 
kinematic endpoints [1,7,12] 

{mf^n^ = ml Rcd (1 - Rbc) (1 - Rab); (1.2) 
(m-,f )' = ml{l- Rcd) (1 - Rbc) ; (1.3) 



( 



{m]}r)' 



mf/^'^y = ml{l-RcD){l-RAB); (1.4) 

' ml{l - i?cD)(l - Rac), for Rcd < Rac, case (1, -), 
ml{l - Rbc){^ - RabRcd), for Rbc < RabRcd, case (2, -), 

(1-5) 

ml{l - i?AB)(l - Rbd), for Rab < Rbd, case (3, -), 

^ ml (l — V Rad ) ^ ) otherwise, case (4, — ). 

Here and below we follow the notation and conventions of Ref. [47], i.e. we write all results 
in terms of an overall mass scale (given by the mass mo oi the heaviest BSM particle D) and 
three dimensionless squared mass ratios 

2 

TfT. 

Rij = ^, i,j e{A,B,C,D}. (1.6) 

Note that there are only three independent ratios in (1.6). We shall take those to be Rab, 
RbCj and Rcd (see Fig. 1), and their definition domain will be the interval (0, 1).^ 

In spite of their transparent theoretical meaning, the set of four endpoints (1.2-1.5) by 
themselves have (justifiably) never been used as the sole basis for a SUSY mass determination 
analysis. This is due to three generic problems, which are all very well known, and are 
separately reviewed in the next three subsections 1.1, 1.2 and 1.3. Our new approach to 
resolving these three problems, and the outline of the rest of the paper are presented in 
Sec. 1.4. 

1.1 Near- far lepton ambiguity 

The first problem is that one cannot diflFerentiate between the "near" and "far" leptons £n 
and if on an event-by-event basis. Since all decays in Fig. 1 are prompt, both leptons point 



^We shall see below that the formulas simplify considerably if we use invariant masses squared instead. 
This distinction is not central to our analysis. 

^As seen in eq. (1.5), at times we shall also utilize one or more of the other three ratios, Rac, Rad and 
Rbd, whenever this will lead to a simplification of the formulas. Of course, the latter three ratios are related 
to our preferred set {Rab, Rbc , Rcd} due to the transitivity property RijRjk = Rik- 
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back to the primary interaction vertex and there is no way to tell which came first and 
which came second. Consequently, one cannot separately construct the individual mj£^ and 
rriji^ invariant mass distributions, whose upper endpoints would be given by (1.3) and (1.4). 
This problem has motivated most of the previous invariant mass studies in the literature, 
beginning with [7], to introduce an alternative definition of the two j£ distributions, simply 
by ordering the two mj^ entries in each event by invariant mass as follows 



mji(io) = mm {mji^,mjef} , 



(1.7) 
(1.8) 



Both of the newly defined quantities mj£(^io) and mj£(^fii) ^'Iso exhibit upper kinematic endpoints 

(m"^?j^-j and Tn^(ni,j\, correspondingly). Since the individual mje(iQ-^ and 



distributions 



are observable, their endpoints are experimentally measurable and can be related to the 
underlying SUSY mass spectrum as follows [7, 12] 



( max y 



m 



'jin I 



m 



[1.9) 



m 



Meq)) 



m 



m 



for (2 - Rab)~'^ < Rbc < 1, case (-, 1), 
for Rab < Rbc < (2 - Rab)~^, case (-,2), 
for < Rbc < Rab, case (—,3); 

, for (2 - Rab)'^ < Rbc < 1, case (-, 1), 

, for Rab < Rbc < (2 - Rab)~\ case (-, 2), (1.10) 



m 



■j4 



\2 

M , for < Rbc < Rab, 



case (—,3); 



where 



K'm)' = ^li^- Rcd) (1 - Rab) (2 - Rab)'' 



(1.11) 



and m^"^ and m^^^ were already defined in (1.3) and (1.4), correspondingly. With this 
approach, the original set of 4 endpoints in eqs. (1.2-1.5) is replaced by 



max max max max 



(1.12) 



In contrast to this conventional approach in the literature, we shall adopt a very different 
attitude towards resolving the problem of the near-far lepton ambiguity. We will do the 
simplest possible thing, namely, we shall do nothing. We shall never ask the question "which 
lepton was in and which one was £fT\ We shall also not use the ordering (1.7,1.8). Instead, 
we shall simply take the two iriji entries in each event, and always treat them in a symmetric 
fashion. For example, any observable invariant mass distribution that we will build out of 
the two measured quantities ruji^ and ruji^ should be invariant under the symmetry 



(1.13) 
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The advantages of our approach may not be immediately obvious at this point, but will 
become clear in the process of our mass determination analysis in Sec. 3 below. 



1.2 Insufficient number of measurements. 

The second problem associated with the original set of four measurements (1.2-1.5), as well 
as the alternative set (1.12), is that the measured endpoints may not all be independent from 
each other. Indeed, there are certain regions of parameter space where one finds the following 
correlation [12] 

In this case, the four measurements (1.12) are clearly insufficient to pin down all four indepen- 
dent input parameters niu, nic, ms and niA- Therefore, one has to measure an additional 
independent endpoint. To this end, it has been suggested to consider the constrained distri- 
bution which exhibits a useful /ower kinematic endpoint "^^^(W-) ['''] 



{m]f;i,^.^y = ^mll (1 - Rab){1 - Rbc){1 + Rcd) (1-15) 



+ 2 (1 - Rac){1 - Rcd) - (1 - Rcd)\/{} + Rab?{1 + Rbc? - 16i?Ac|- 

The distribution mj££(^Qy^j is nothing but the usual rriju distribution over a subset of the 
original events, subject to the additional dilepton mass constraint 



^max 



V2 



< mu < mf^"-'' . (1.16) 



In the rest frame of particle B, this cut implies the following restriction on the opening angle 
9 between the two leptons [6] 

e>^, (1.17) 

thus justifying the notation for mji£(^Qy^y 

The advantage of the "threshold" endpoint measurement (1.15) is that it is always in- 
dependent of the other four measurements in (1.12). As a result, it would appear that the 
enlarged set of five kinematic endpoint measurements 

^max ^max ^max ^max ^min ( -i -io\ 

^et ^^ju '^je{io)'^je{hi)^^jee{e>^) \^-^°) 

should be in principle sufficient to determine all four unknown masses (see, however, Ref. [47] 
and Section 1.3 below). 

Unfortunately, the "threshold" (1.15) also suffers from certain disadvantages, which are 
mostly of experimental nature. It is generally expected that the experimental precision on 
the determination of the lower kinematic endpoint (1.15) will be rather inferior compared 
to the precision on the other four upper kinematic endpoints (1.12) [12]. There are several 
generic reasons for such a pessimistic attitude. First, the region in the 'mju{e>^) distribution 
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near its lower endpoint (1.15) is rather sparsely populated, resulting in a shallow edge and 
sizable statistical errors. To make matters worse, the niji^^gy^-^ distribution near its lower 
edge is a convex function [19], which makes it even more difficult to tell where the signal 
ends and the tails from various sources begin [12]. Finally, the low mass region of almost any 
invariant mass distribution in SUSY is generally associated with larger SM (as well as SUSY 
combinatorial) backgrounds compared to its high mass counterpart. 

Overall we find all these disadvantages sufficiently convincing so that we will drop the 
measurement (1.15) altogether and will never use it in the course of our analysis in Sec. 3 
below. We will be justified in doing so, since the linear dependence problem (1.14), which 



has plagued previous studies and was the prime motivation for introducing the "i^^^^^^i) 
measurement in the first place, will have no effect on our analysis. In fact, we will not be using 
the endpoint measurement JT^^^^j) (for the reasons given in the previous subsection 1.1) and 
we will not be using the endpoint measurement rn!^^^ (for the reasons given in the following 
subsection 1.3). Once these two problematic measurements are removed from consideration, 
the linear dependence problem (1.14) does not arise, and the "threshold" measurement (1.15) 
is not central to the analysis any more. 

1.3 Parameter space region ambiguity 

The third problem with the conventional set of measurements (1.18) is immediately obvious 
from the defining equations (1.5), (1.9) and (1.10) for the kinematic endpoints m^"^, mnJ^^i^y 
and iT^Ye(hi)^ correspondingly. One can sec that the relevant expressions are piecewise-defined 
functions, i.e. they depend on the values of the independent variables uia, mB, mc and uid- 
For example, there are four different cases for m^^^ , and three different cases for the pair of 
("^^fio)' "^^(/li))- Altogether, these give rise to 9 different cases^ which must be separately con- 
sidered [12,47]. Of course, this represents a problem, since the masses are a priori unknown, 
and it is not clear which case is the relevant one. Barring any model-dependent assumptions, 
one is forced to consider all possibilities, obtain a solution for the spectrum, and only at the 
very end, test whether the solution falls within the parameter space applicable for the case at 
hand. This procedure may often result in several alternative solutions [12,47,78-82]. In fact, 
Ref. [47] recently proved that there exists a sizable parameter space region in which even the 
full set of measurements (1.18) would always yield two alternative solutions, even under ideal 
experimental conditions. The problem is further exacerbated by the inevitable experimental 
errors on the measurements (1.18), which would allow for an even larger number of "fake" or 
"duplicate" solutions [47,79,80]. 

Having identified the root of the duplication problem as the piecewise definition of the 
mathematical formulas in (1.5,1.9,1.10), our solution to the problem will be again very sim- 
ple and conservative. We will simply avoid using any kinematic endpoints which are given 
in terms of piecewise-defined expressions. This requirement automatically eliminates from 
consideration the three conventional endpoints mJ^\^^ ■, ^Yi{lo)^ ''^jiihi)' ^^^^^ already 
gave up on mYule^s.^ ™ the previous subsection, this leaves m^"^ as the only measurement 



''The remaining 3 cases are always unphysical [12]. 
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out of the conventional set (1-18) that we shall use in our analysis. This is perhaps the most 
drastic difference between our approach and all previous studies in the literature. 

1.4 Posing the problem 

In the previous three subsections we discussed each of the three generic theoretical^ problems 
with the previous applications of the kinematic endpoint method for mass determination. 
We arc now ready to explicitly formulate our main goal in this paper. We aim to design a 
method for measuring the masses of the particles in the decay chain of Fig. 1, which is based 
on kinematic endpoint information, and satisfies the following requirements: 

• It does not make use of any kinematic endpoints whose interpretation is ambiguous, 
i.e. whose expressions in terms of the physical masses are piecewise-defined functions. 

• It does not make use of any lower kinematic endpoints such as the "threshold" 'm-Y£i{e>^)^ 
due to the experimental challenges with such measurements. 

• It relies solely on 1-dimcnsional distributions, unlike the methods recently advertised 
in [18,46,47,81,82], which utilize 2-dimensional correlation plots. While the latter 
do provide a wealth of valuable information, they also typically require more data in 
order to obtain good enough statistics for drawing any robust conclusions from them. 
In contrast, the one-dimensional distributions should be available rather early on, and 
with sufficient statistics for endpoint measurements. 

As already alluded to in the previous subsections, the first two requirements already eliminate 
four out of the five conventional inputs (1.18). Obviously, we will need to find a way to replace 
those with an alternative set of kinematic endpoint measurements which nevertheless satisfy 
the above requirements. In Section 2 we introduce and investigate a new set of invariant 
mass variables whose upper endpoints can be useful for our analysis. Then in Section 3 we 
outline our basic method, which makes use of some of these new variables. We illustrate our 
discussion in Section 4 with two numerical examples: the LMl and LM6 CMS study points. 
Section 5 is reserved for our conclusions. In Appendix A we supply the analytic expressions 
for the shapes of the 1-dimensional invariant mass distributions used in our main analysis 
in Sec. 3.1. Those results can be useful in improving the precision on the extraction of the 
kinematical endpoints. 

*In addition, there are problems which are of experimental nature, e.g. identifying the correct jet and 

the correct lepton pair resulting from the decay chain in Fig. 1. There exists a set of standard experimental 
techniques which are aimed at overcoming these problems, e.g. the opposite flavor subtraction for the two 
leptons and the mixed event subtraction for the jet [83]. Wrong U and j£ pairings can also be identified 
and a posteriori removed whenever an invariant mass entry for rnei, rriji or rnjtt exceeds the corresponding 
kinematic endpoint m^(°'^, m'Jl'f^j^i'^ or m'Jll^. In what follows we shall assume that those preliminary steps 
have already been done and the samples we are dealing with have already been appropriately subtracted to 
remove the combinatorial background. 
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2. New vairiables 



In this section we propose a new set of invariant mass (squared) variables. As already 
explained in the Introduction, our variables should be composed of rnj^^ and m'jg^ in a 
symmetric way, in accordance with (1.13). Consequently, any plotting manipulations or 
mathematical operations involving m^^^ and friji^ should obey the symmetry implied by 
eq. (1.13). 

2.1 The union rrij^^ U ruj^^ 

We begin with the simplest case, where we postpone applying any mathematical operations 
to iTij^^ and 'mje^j and instead simply plot them. The requirement of eq. (1.13) implies that 
the only possibility is to place both of them together on the same plot, in essence forming 
the union 

of the individual fn'j^^ and fn'jij, distributions. Since each individual distribution is smooth 
and has a kinematic endpoint, the same two kinematic endpoints should be visible on the 
combined distribution "rn^^u^j as well^. We shall denote the larger of the two endpoints with 



and the smaller of the two endpoints with 



MI^*,?.^! ^ = max I (m^;^^) ^ , (m^^^) ^ | (2.2) 



The newly introduced quantities MV^^ and 'ti^^'^) are nothing but the usual kinematic end- 
points m^"^ and m'j}'^^, given by (1.3) and (1.4), correspondingly. Of course, at this point 
we do not know which is which, and we have an apparent two-fold ambiguity: we can have 
either 

3l{u) - ^jin ' ^jl{u) - ^jif ' " > HbC, (2.4j 

or 

jl(u) - ^jif ' ^jl{u) - ^jin ' " ^AB < J^BC- (^-5) 

Notice that both (2.2) and (2.3) are officially upper kinematic endpoints, and thus satisfy our 
basic requirements. 

The benefits of our alternative treatment (2.1) in response to the near-far lepton am- 
biguity problem of Sec. 1.1, arc now starting to emerge. With the conventional ordering 
(1.7,1.8) one has to deal with a three-fold ambiguity in the interpretation of the endpoints 
''^Jello) ^j^ffet)' ^ seen in eqs. (1.9,1.10). Instead, the simple union (2.1) leads only to 
the two-fold ambiguity of eqs. (2.4,2.5). More importantly, the analysis of Sec. 3.1 below will 
reveal that in spite of the remaining two-fold ambiguity in eqs. (2.4,2.5), one can nevertheless 
uniquely determine all three of the masses mo, rnc and m^! We consider this to be one of 
the important results of this paper. 

®For specific numerical examples, refer to Sec. 4. 
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2.2 The product mji^ x rriji^ 

In the remainder of this section, we shah construct new invariant mass squared variables out 
of the two entries rn^f^ and 'f^'ji^.^ simply by applying various mathematical operations on 
them in a symmetric fashion. We begin with the product 

^%(P) = '^jtn'mjef (2.6) 

whose endpoint is given by 

' ^ml{l- RcdWI - Rab, for Rbc < 0.5, 



2 

niT 



(2.7) 



.^(1 - Rcd)VRbc{1 - Rbc){1 - Rab), for i?Bc > 0.5. 
Unfortunately, this endpoint also turns out to be piecewise-defined, thus failing one of our 
basic requirements from the Introduction. Therefore we shall not use this endpoint in the 
course of our analysis. 

2.3 The sums mjl^ + m^"^ 

Another possibility is to consider various sums, for example w-j^^ + "^f^^ i^jl-n + '^i^/)^) 
as originally proposed in [18]. Here we generalize the discussion in [18] and introduce a whole 
set of new variables, m'^^^^^{a), labelled by the continuous parameter a, which are defined as 

m%is) («) = + m^if) " ■ (2-8) 

Since a is a continuous parameter, in principle there are infinitely many mj£(^g) variables! 
Notice that the conventional variables '^^^^(jo) and ^rz^^^^j^ from (1.7) and (1.8) are also included 
in our set, and are simply given by 

"^jV) = "i?€(s)(-oo)' (2-9) 

^"jeihi) = '^'je{s){'^) ■ (2-10) 

We see that our new set (2.8) is a very broad generalization of the conventional definitions 
(1.7) and (1.8), which just correspond to the two extreme cases a = ±00. Of course, the user 
is free to choose a at will, and any finite value of a will lead to a new variable i^^£(^s) ('^) ■ 

In order to make the new variables m'^^^^^{a) useful for mass spectrum studies, we need 
to provide the formulas for their kinematic endpoints (w-^('f)(a))^. These formulas are easy 
to derive, using the results from [47], and we present them in the next two subsections, where 
it is convenient to consider separately the following two cases: a > 1 (in Sec. 2.3.1) and 
a < 1, but a 7^ (in Sec. 2.3.2). 

2.3.1 Kinematic endpoints of 'm'^£(^s^{a) with a> 1 

When one chooses a value of a > 1, the "^|^(5)(ck) endpoint is given by the following expression 



V 



Rab < I- {I- Rbc) (1 - 



h"^rr)(«>i)) - ^ . ^ (2.11) 

^ (m^T'^(a)j , Rab >l-il- Rbc) (1 - R%c) 
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where m^"^ was already defined in (1.4), and m 
quantity 



max 



(a) is a newly defined, a-dependent 



{mfr{a)f = ml{l - Rcd) R%ci^ - RabT + (1 - RbcT 



(2.12) 



As a cross-check, one can verify that in the limit a — oo the expression (2.11) reduces to 
(1.10), in agreement with (2.10). In that case, the upper line in (2.11) corresponds to options 



(-, 1) and (-,2) in (1.10), where m™^^.) 
option (—,3) in (1.10), where m^"^^^ — 



m 



■3h 



m 



■jin 



' , while the lower line in (2.11) corresponds to 
Unfortunately, just like the product endpoint 
(2.7), the endpoint (2.11) is in general piecewise-defined, and does not meet our criteria. 

However, there is one important exception, namely the case of a = 1, in which we do get 
a singly defined function. According to the general definition (2.8), m'^^^^^{a = 1) is simply 
the sum of the two entries in each event: 



[a 



Using the identity 



(2.13) can be equivalently rewritten as 



i£(.)(« = 1) = • 



(2.13) 



(2.14) 



(2.15) 



To find the expression for its endpoint, one can set a = 1 in (2.11), and then realize that the 
logical condition for executing the upper line becomes Rab 0, which is impossible, since the 
mass ratios Rij in (1.6) are always positive definite. Therefore, the endpoint m^^^-j(a = 1) is 
always calculated according to the lower line in (2.11), which results in [18] 



(2.16) 



Note that this endpoint is perfect for our purposes since the formula (2.16) is always unique, 
i.e. it is independent of the parameter space region. The variable m|^^^^(a = 1) will thus 
play a crucial role in our analysis below. 

2.3.2 Kinematic endpoints of "^|£(5)(a) with a < 1 and a 7^ 

Finally, in the case when cc < 1, but a 7^ 0, the f^^m^g-^ip) endpoint is given by the following 
expression 



' (m-/-(a))', 
m\{,\-RcD) \\^{\-Rab)^- 



Rbc > 



, Rbc < 



i + {i-RABy 



1 + {1-Rab)~ 

(2.17) 

where m^^^{a) was already defined in (2.12). Again as a cross-check, one can verify that 
in the limit a — —00 the expression (2.17) reduces to (1.9), in agreement with (2.9). In 
the a — —00 case, the upper line in (2.17) corresponds to option (—,1) in (1.9), where 
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"^^(io) ~ "^j^"^' while the lower line in (2.17) corresponds to options (—,2) and (—,3) in 
(1.9), where fn^'^j^i^-^ = ^"i{eq)' U^ifo'^tunatcly, the cndpoint function (2.17) is again piecewise- 
defined, and does not meet one of our basic criteria spelled out in the introduction. 

In passing, we note that the special case of a = | , which involves the linear sum of the 
two masses ^ 

was previously explored in [18,84]. In that case, from (2.17) we find for its endpoint 



mlil - Rcd) [^/RBc{l - Rab) + VI - RBcf , Rbc > ^ 



-Rab 
2-Rab ' 



mlil - Rcd){2 - Rab), Rbc < ^eIH • 

(2.19) 



2.4 The difference {rnj^^ — rn^^^ \ 

Finally, one can also consider a set of variables which involve the absolute value of differences 
between and nT^j^^. In analogy with (2.8), we can define another infinite set of variables 



(2.20) 



Once again, the user is free to consider arbitrary values of a. However, this freedom is 
redundant, when it comes to the issue of the kinematic endpoints of the variables in (2.20). 



It is not difficult to see that the endpoints of rn^^f/j) (a) are always given by 



(2.21) 

and are in fact independent of a! Therefore, for the purposes of our discussion, it is sufficient 
to consider just one particular value of a. In the following we shall only use a = 1: 



"^jm^oi = 1) = ruj^^ - rriji^ , (2.22) 

which is the analogue of 'm^£(^g^{c( = 1) defined in (2.13). 

The result (2.21) implies that the endpoint of (2.22) does not contain any new amount of 
information, which was not already present in the two kinematic endpoints MJ^y^ and rnj^i^^^ 
discussed in Sec. 2.1. Nevertheless, the independent measurement of (m™|^(l))^ can still be 
very useful, since it will mark the location of (MJJ'^'^p^ on the fn'^^^^ distribution. Then one 
will be looking for the second endpoint (mj^^^-j)'^ to the left, i.e. in the region of smaller fn^K^^-j 
values. 

This completes our discussion of the new invariant mass variables and their kinematic 
endpoints. For our basic proof-of-principle measurement technique presented in the next 
Section 3.1, we shall use only three of them, namely M"^'^^y ''^^je{u)^ "''j^(s)('-* ~ ^)' 
However, the remaining variables are in principle just as good, their only disadvantage being 
that they failed our arbitrarily imposed condition at the beginning that the endpoint functions 
should all be region independent. Of course, one could, and in fact should, use all of the 
available kinematic endpoint information, which in a global fit analysis can only increase the 
experimental precision of the sparticle mass determination. 
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3. Theoretical analysis 



3.1 Our method and the solution for the mass spectrum 

Our starting point is the set of four measurements 

max lijmax max rnmax ( _ -,\ 

in place of the conventional set (1.18). It is easy to verify that the measurements (3.1) are 
always independent of each other, and thus never suffer from the linear dependence problem 
discussed in Section 1.2. 

Given the set of four measurements (3.1), it is easy to solve for the mass spectrum. To 
simplify the notation, we introduce the following shorthand notation for the endpoints of the 
mass squared distributions 

L^{mTen\ M^{mJII:^\ m^(mf^:^\ S ^ [mf,^^{a = l)f (3.2) 
The solution for the mass spectrum is then given by 

2 _ Mw{L + M + rn-S) _ 

2 MmL 

™^ = {M + m-Sf ^'-'^ 



(M+m-sy^ ' if Rab > Rbc, 

mL(S—m) -c TD ^ o 

(M+m-s)^ ' if Rab < Rbc; 



(3.5) 



^ _ L{S-m){S-M ) 
{M + m- Sf 

It is easy to verify that the right-hand side expressions in these equations are always positive 
definite, so that one can safely take the square root and compute the linear masses mu, mc, 
ruB and uia- Notice that in spite of the two-fold ambiguity (2.4,2.5), the solution for m/j, 
mc and mA is unique! Indeed, the expressions for mn, mc and are symmetric under 
the interchange M m. The remaining two-fold ambiguity for is precisely the result 
of the ambiguous interpretation (2.4,2.5) of the two "^^^(u) endpoints, and is related to the 
symmetry under (1.13), or equivalently, under the interchange 

Rab ^ Rbc ■ (3.7) 

In the next subsection we discuss several ways in which one can lift the remaining two-fold 
degeneracy for ms which is due to (3.7). 

Notice the great simplicity of this method. The expressions for (3.3), (3.4) and (3.6) are 
region independent and therefore one does not have to go through the standard trial and 
error procedure involving the 9 parameter space regions {Nju,Nj£) [12,47] associated with 
the various interpretations of the endpoints m^^^ , "^^(io) "^^(^«)- 
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Figure 2: Comparison of the predictions for the kinematic endpoints rn™g^'^-^{a) of the real and fake 
solutions, as a function of 4> = arctana (in units of t:), for the two examples discussed in detail in 
Section 4: (a) the LMl CMS study point and (b) the LM6 CMS study point. In each panel, the 
prediction of the real (fake) solution is plotted in red (blue). The vertical dotted line indicates the 
case of (j) = ^ (a — 1), for which the two solutions give an identical answer, marked with a green dot. 
The horizontal dotted lines show the corresponding asymptotic values m^"^--) and m^^j^^y obtained 
at a ^ ±oo ((j) ±f j- 



3.2 Disambiguation of the two solutions for hib 

The method outhned in Sec. 3.1 allowed us to find the true masses of particles A, C and 
D, but yields two separate possible solutions for the mass of particle B. We shall now 
discuss several ways of lifting the remaining two- fold degeneracy for niB- 

3.2.1 Invariant mass endpoint method 

One possibility is to use an additional measurement of an invariant mass endpoint. Indeed, as 
shown in Sees. 1 and 2, there are still quite a few one-dimensional invariant mass distributions 
at our disposal, which we have not used so far. Those include the conventional distributions 
of m|^^, and m^jn^^iy ^^^^ "^^"^ distributions m|^(p), and m^^^^^ 

(1) 

which we introduced in Sec. 2. Which of them can be used for our purposes? Note that the 
duplication in (3.5) arose due to the symmetry (3.7), so that any kinematic endpoint which 
violates this symmetry will be able to distinguish between the two solutions. 

Let us begin with the conventional distributions "t-^^^, ^'^^(/o); ^^{hi) ^'jU{e>^)^ 
whose endpoints we did not use in our analysis so far. It is easy to check that m^"^, iT^Ye{hi) 
and ^Yeel^ey-) invariant under the interchange (3.7) and cannot be used for discrimination. 
However, m"^^^-^ is not symmetric under (3.7) and can do the job. In fact, one can show that 
the two duplicate solutions for m b always give different predictions for m^^f^-^ . 

More importantly, many of our new variables from Sec. 2 can provide an independent 
cross-check on the correct choice for the solution. For example, the kinematic endpoint (2.7) 
of the product variable also violates the symmetry (3.7) and distinguishes among the 

^The only exception is the trivial case of Rab = Rbc, but then the two solutions for ms coincide, and 
ms is again uniquely determined. 
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two solutions. The infinite set of variables '^^^(s)(Q!) can also be used, and for almost the whole 
range of a < 1. To see this, in Fig. 2 we compare the predictions for the kinematic cndpoints 
m^"^^(a) of the real and fake solutions, for the two examples discussed in detail in Section 4: 
(a) the LMl CMS study point and (b) the LM6 CMS study point. The corresponding mass 
spectra are listed in Table 1 below. For convenience, we plot versus the parameter 

(f) = arctana , (3-8) 

which allows us to map the whole definition domain (—00,00) for a into the finite region 
(— ^, ^) for (/). Fig. 2 shows that for most of the allowed (j) range, the two solutions predict 
different values for the kinematic endpoints rn^^^s^{a). In fact, for < ^, the two predictions 
are always different, apart from the trivial case of = (a = 0). Even for ^ > ^, there 
still exists a range of ^, for which, at least theoretically, a discrimination can be made. The 
predictions are guaranteed to coincide only for = | (a = 1) (as they should, see (3.1)), and 
for a certain range of the largest possible values of 0. 

3.2.2 Invariant mass correlations 

Another way to resolve the twofold ambiguity in our solution (3.5) is to simply go back to 
the original measurements of M^^^^ and rri^^^^ and already at that point try to decide which 
of the two measured rUjK^u) endpoints is m^"^ and which one is m^"^. As already discussed 
in [18,46], this identification is in principle possible, if one considers the correlations which are 
present in the two-dimensional distribution versus mfi. The basic idea is illustrated 

in Fig. 3, where we show scatter plots of mj^(„) versus m^^, for the two examples used in 
Fig. 2 and discussed in detail later in Section 4. Fig. 3(a) (Fig. 3(b)) shows the result for the 
real (fake) solution corresponding to the LMl study point, while Figs. 3(c) and 3(d) show 
the analogous results for the LM6 study point. In each plot we used 10,000 entries, which 
roughly corresponds to 20 fb"^ (200 fb"^) of data for the actual LMl (LM6) SUSY study 
point. Here and below we show the ideal case where we neglect smearing effects due to the 
finite detector resolution, finite particle widths and combinatorial backgrounds. All of our 
plots are at the parton level (using our own Monte-Carlo phase space generator) and without 
any cuts. Notice that in order to avoid dealing with the large numerical values of the squared 
masses, we use a quadratic power scale on both axes, which allows us to preserve the simple 
shapes of the scatter plots when plotting versus the linear masses themselves. 

Fig. 3 shows that the combined distribution is simply composed of the two separate 

distributions rn^jf^^ and "t-^^^ , but they are correlated differently with the dilepton distribution 
m^^. In particular, let us concentrate on the conditional maxima mJJgi^ {mu) and ^^"^(m^^), 
i.e. the maximum allowed values of rrij^^ and m^^^, respectively, for a given fixed value 
of mu [18,46]. A close inspection of Fig. 3 shows that the values of m^^^ and m^^ are 
uncorrelated, and as a result, the conditional maximum {mu) does not depend on mu. 

In turn, this implies that the endpoint value {m^^)"^ given in (1.3) can be obtained for any 
m%: 

n = [mfCf = [mfCimuf = ml{l - Rcd) (1 - Rbc), Vm« G [0,m??-] . (3.9) 
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Figure 3: Predicted scatter plots of mj^^^-^ versus mu, for the case of the real and fake solutions 
for each of the two study points LMl and LM6: (a) the real solution LMl; (b) the fake solution 
LMl'; (c) the real solution LM6; and (d) the fake solution LM6'. The red solid horizontal (blue 
dashed inclined) line indicates the conditional maximum m^""^ (mu) (m^'^^{mu)) given by eq. (3.9) 
(eq. (3.10)). Each panel contains 10,000 entries. The results shown here are idealized in the sense that 
we neglect smearing effects due to the finite detector resolution, finite particle widths and combinatorial 
backgrounds. Notice the use of quadratic power scale on the two axes, which preserves the simple 
shapes of the scatter plots, even when plotted versus the linear masses mj^u) and mu. 



Because of (3.9), the shape of the m-^^^ versus mj^ scatter plot is a simple rectangle [18,46]. 
This is confirmed by the plots in Fig. 3, where the (red) horizontal solid line indicates the 
constant value (3.9) for the conditional maximum ra'^^{mu). 

In contrast, the values of f^^jij and are correlated. The conditional maximum 
m^^^ {mit) does depend on the value of m^^ as follows: 

(m5^7(m«))' =p+ ^m|, (3.10) 
where we introduce the shorthand notation used in [47] 

/ = {rnjt;)^ = ml{l- Rcd) (1 - Rab), (3.11) 

P = RBcf = ml{l- Rcd) Rbc (1 - Rab). (3.12) 



- 15 - 



The absolute maximum of n^'jij:^ which is given by (1.4) and denoted here by /, can only be 
obtained when mj^ itself is at a maximum [18,46]: 

(3.13) 

On the other hand, the conditional maximum rnJ^^irnu) obtains its minimum value at 
•mff^ = and corresponds to [18, 46] 

'</. (3.14) 

Eqs. (3.13,3.14) imply that the shape of the "^^^^ versus m^^ scatter plot is a right-angle 
trapezoid. This is confirmed by the plots in Fig. 3, where we mark with a (blue) dashed line 
the conditional maximum (3.10). With sufficient statistics, this difference in the kinematic 
boundaries may be observable, and would reveal the identity of rnJ^^ and m^"^ [18,46]. Once 
the individual m^"^ and m^"^ are known, the solution for the mass spectrum is unique - see 
e.g. Appendix A in [47]. Of course, in cases where p ~ /, namely i?BC' ~ 1, it may be difficult 
in practice to tell which of the two boundaries in the scatter plot is inclined and which one 
is horizontal^. One example of this sort is offered by point LM6, which has Rbc = 0.91 and 
leads to a rather fiat rn^-f^^{mie) function, as seen in Fig. 3(c). 

An alternative and somewhat related method will be to investigate the shapes of the 
one-dimensional distributions themselves [85]. In Appendix A we provide the analytical 
expressions for the shapes of the four invariant mass distributions m|^, m?^^^^y m^^^^j(l) and 
m|^^^-j(l) used in our basic analysis from Sec. 3.1. Given what we have already seen in 
Fig. 3, it is not surprising that the true and the fake solutions predict different shapes for 
the one-dimensional distributions as well. In the LMl and LM6 examples considered below 
in Sec. 4, this difference is particularly noticeable for the "T'j£((i)(l) distributions 

(see Figs. 4(b), 4(d), 5(b) and 5(d)), and can be tested experimentally. 

3.2.3 Mt2 endpoint method 

Let us note that if we identify particle A with the LSP, we have a rather peculiar situation, in 
which we know the LSP mass niA-, and we are unsure about the NLSP mass Mb-, for which we 
have to choose among two alternatives. This goes against the common lore which considers 
the LSP mass (in this case m^) to be the least constrained among the masses appearing in 
the decay chain in Fig. 1. For example, the method of the Cambridge Mt2 variable [2,8] 
treats the LSP mass as a continuous unknown parameter. At this point of our analysis we 
already know the LSP mass, and we can use this knowledge to our advantage. For example, 
if we can collect a sufficient number of events of B pair-production, we can apply the idea 
of Mt2 for the B ^ A decay as in the original Mt2 proposal [2]. When we use for the 
trial LSP mass the known true value of ruA given by (3.6), the kinematic endpoint of the 
Mt2 distribution will reveal the correct value of the mass uib oi the parent particle B, thus 
selecting the true solution in (3.5). 

*A separate problem, which arises in the case of p ~ /, will be discussed below in Sec. 4.1. 
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As emphasized in Ref. [43], the Mt2 endpoint method does not necessarily rely on A 
being the LSP (i.e. the very last partiele in the decay chain) or B being the "grandparent" 
(i.e. the very first particle in the decay chain). For example, suppose that A decays further. 
In that case, one simply needs to apply the more general "subsystem" variable Mj^^^''^^ [43] 
with A being the "child" particle: c = A. Similarly, the two B particles do not have to be the 
two grandparents initiating the decay chains: it is sufficient to consider m!^^^''^^ with p = B 
and arbitrary n [43]. Finally, for the purposes of selecting the correct solution in (3.5) it is 
also possible to apply the subsystem variable M^f^^''^^ in a different way, where B is the child, 
and the parent is either D or C. In this case, we know the parent mass, which is respectively 
given by (3.3) or (3.4), and we are asking the question, which of the two test masses in (3.5) 
gives the correct answer for the Mt2 endpoint. 



4. Numerical examples 

We shall now illustrate the ideas of the previous section with two specific numerical examples: 
the LMl and LM6 SUSY study points in CMS [83]. The mass spectra at LMl and LM6 are 
listed in Table 1. Point LMl is similar to benchmark point A (A') in Ref. [86] (Ref. [87]) 
and to benchmark point SPSla in Ref. [88]. Point LM6 is similar to benchmark point C (C) 
in Ref. [86] (Ref. [87]). The table also lists the corresponding duplicate solutions LMl' and 
LM6', which are obtained by interchanging Rbc ^ Rab, or equivalently, by replacing the 
mass of B via 

/ rnAUic , . 

mB^ms = . (4.1) 

rriB 

It is interesting to note that LMl and LM6 represent both sides of the ambiguity (3.7): 
at LMl, we have Rab > Rbc and correspondingly, m^^^ > m^"^ and (2.4) applies. On the 
other hand, at LM6 we have Rab < Rbc and m^"^ < m^"^, so that (2.5) applies. Another 
interesting difference is that at LMl particle B is the right-handed slepton Ir, while at LM6 
the role of particle B is played^ by the left-handed slepton II. Of course, to the extent that 
we are interested in kinematical features, this difference is not relevant, and particle B of the 
LM6 spectrum may very well have been the right-handed slepton instead. 

4.1 Mass measurements at points LMl and LM6 

Given the mass spectra in Table 1, it is straightforward to construct and investigate the 
relevant invariant mass distributions. For the purposes of illustration, we shall ignore spin 
correlations, referring the readers interested in those effects to Refs. [52,58,71]. We are 
justified to do so for several reasons. First, our method relies on the measurement of kinematic 
endpoints, whose location is unaffected by the presence of spin correlations. Second, in 
the case of supersymmetry (which is really what we have in mind here), particle S is a 
scalar, which automatically washes out any spin effects in the mj^ and m^^^ distributions. 
Furthermore, if particles D and their antiparticles D are produced in equal numbers, as would 



^Although the right-handed slepton £r is also kinematically accessible at point LM6, the wino-like neu- 
tralino £2 decays much more often to £l as opposed to £r. 



-17- 



Variable 


LMl 


LMl' 


LM6 


LM6' 


niA (GeV) 


94.9 


158.15 


mB (GeV) 


118.9 


143.35 


291.0 


165.65 


mc (GeV) 


179.6 


304.8 


ran (GeV) 


561.6 


861.9 


Rab 


0.6370 


0.4383 


0.2954 


0.9115 


Rbc 


0.4383 


0.6370 


0.9115 


0.2954 


RCD 


0.1023 


0.1251 


mf^°^ (GeV) 


81.10 


76.12 


(GeV) 


398.8 


676.8 


(GeV) 


320.6 


239.8 


"^"£"f)(« = 1) (GeV) 


451.8 


689.2 


m™^^ (GeV) 


451.8 


689.2 


^mo>^) (GeV) 


215.2 


176.4 


"^,"^(.".) (GeV) 


398.8 


676.8 


m™^f)(a = 2) (GeV) 


406.6 


398.8 


676.8 


677.0 


m™^f)(a = 1.5) (GcV) 


417.9 


402.5 


676.8 


678.4 


m™^f)(a = 0.5) (GeV) 


611.0 


638.9 


886.0 


807.1 


mf^'(;s){(x = -0.5) (GeV) 


142.9 


159.7 


174.9 


138.0 


^Te(s)(^ = -1) (GeV) 


200.1 


225.9 


224.8 


184.8 


"^J^(Zo) (GeV) 


274.6 


319.1 


239.8 


229.9 


m^^^) (GeV) 


292.0 


319.4 


393.7 


310.9 


m^'fT (GeV) 


398.8 


320.6 


239.8 


676.8 


^marr (GeV) 


320.6 


398.8 


676.8 


239.8 



Table 1: The relevant part of the SUSY mass spectrum for the LMl and LM6 study points. The 
corresponding duplicated solutions LMl' and LM6' are obtained by interchanging Rbc ^ Rab as in 
(3.7). In the table we also list the corresponding values for various invariant mass endpoints. The first 
four of those represent our basic set of measurements (3.1) discussed in detail in Section 4.1, while the 
last two (mj^"^ and mj^"^) are not directly observable. The remaining invariant mass endpoints are 
considered in Section 4.2. In the case of mj]"^) (a), we show several representative values for a. For the 
complete a variation, refer to Fig. 2. Recall that ?Ti^"f)(+oo) = m^^'^^j and to^('s^)(— oo) = rn^^'^i^y 

be the case if the dominant production is from gg and/or qq initial state, any spin correlations 
in the rn^j^^ distribution are also washed out. Under those circumstances, therefore, the pure 
phase space distributions shown here are in fact the correct answer. 

We begin our discussion with the four invariant mass distributions m^^, to^^^^-^, m'^^^^^{a = 
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Figure 4: One-dimensional invariant mass distributions for tire case of LMl (red solid lines) and 
LMl' (blue dotted lines) spectra. The kinematic endpoints (3.1) used in our analysis in Section 3.1 
can be observed from these distributions as follows: m^^"^ is the upper kinematic endpoint of the mu 
distribution in panel (a); MJ^^^-^ is the absolute upper kinematic endpoint seen in both the combined 
mjK^u) distribution in panel (b), or the difference distribution mjg(^d.)i^) in panel (d); m.J'^'jy^ is the 
intermediate kinematic endpoint seen in panel (b); and = 1) is the upper kinematic endpoint 

of the mjHs'j{a = 1) distribution in panel (c). 



1) and fn^K^^-^ict = 1), which form the basis of our method outhned in Sec. 3.1. Fig. 4 (Fig. 5) 
shows those four distributions for the case of study point LMl (LM6). In each panel, the red 
(soHd) histogram corresponds to the nominal spectrum (LMl or LM6), while the blue (dotted) 
histogram corresponds to the "fake" solution (LMl' or LM6'), which is obtained through 
the replacement (4.1). For all figures in this section, we use the same 4 samples of 10,000 
events each, which were already used to make Fig. 3. Notice our somewhat unconventional 
way of filling and then plotting the histograms in this section. First, we show differential 
distributions in the corresponding mass squared, i.e. dN/dm?. This is done in order to 
preserve the connection to the analytical results in Appendix A, which are written the same 
way. More importantly, the shapes of the one-dimensional histograms are much simpler in 
the case of dN jdrr? as opposed to dN/dm [52,58,71]. In the next step, however, we choose to 
plot the thus obtained histogram versus the mass itself rather than the mass squared. This 
allows one to read off immediately the corresponding endpoint and compare directly to the 
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values listed in Table 1. It also keeps the rc-axis range within a manageable range. However, 
since the histograms were binned on a mass squared scale, if we were to use a linear scale on 
the rc-axis, we would get bins with varying size. This would be rather inconvenient and more 
importantly, would distort the nice simple shapes of the dN/dm? distributions. Therefore, we 
use a quadratic scale on the x-axis, which preserves the nice shapes and leads to a constant 
bin size on each plot. 

Figs. 4 and 5 illustrate how each one of the measurements (3.1) can be obtained. For 
example, m^"^ is the classic upper kinematic endpoint of the mu distributions in Figs. 4(a) 
and 5(a). This endpoint is very sharp and should be easily observable. M^^^-^ is the absolute 
upper kinematic endpoint seen in the combined mj^(y^ distribution in Figs. 4(b) and 5(b). 
Notice that the same endpoint can independently also be observed as the absolute upper 
kinematic limit of the difference distributions mji[d){^) shown in Figs. 4(d) and 5(d). The 
fact that there are two independent ways of getting to the endpoint M^^^-^ should allow for 
a reasonable accuracy of its measurement. Upon closer inspection of the combined mj^f^u) 
distribution in Figs. 4(b) and 5(b), we also notice the intermediate kinematic endpoint n-^"^) 
seen around 320 GeV in Fig. 4(b) and around 240 GeV in Fig. 5(b). Finally, m'^^^^^{a = 1) is 
the upper kinematic endpoint of the mj£(^g^{a = 1) distribution shown in Figs. 4(c) and 5(c). 
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It is also rather well defined, and should be well measured in the real data. 

At this point we would like to comment on one potential problem which is not immedi- 
ately obvious, but nevertheless has been encountered in practical applications of the invariant 
mass technique for SUSY mass determinations [85]. It has been noted that in the case of 
p ^ f (see eqs. (3.11,3.12)), the numerical fit for the mass spectrum becomes rather unstable. 
Given our analytical results in Sec. 3.1, we are now able to trace the root of the problem. 
Notice that p ^ f implies that Rbc ~ 1- In this limit, from eqs. (1.2), (1.3), (1.4) and (2.16) 
we find 

lim (L) = 0, lim (ra) = 0, lim {M + m - S) = 0. (4.2) 

iiBC— *1 RbC-^^ RbC-^^ 

This means that the functions (3.3-3.6) giving the solution for the mass spectrum will all 
behave as j^, and, given the statistical fluctuations in an actual analysis, will have very 
poor convergence properties. We note that this problem is not limited to our preferred set 
of measurements (3.1) and is rather generic, but has been missed in most previous studies 
simply because the case of i^^c ~ 1 was rarely considered. 

Figs. 4 and 5 reveal that, as expected, the real (red solid lines) and fake (blue dotted 
lines) solutions always give identical results for our basic set of four cndpoint measurements 
(3.1). This is by design, and in order to discriminate among the real and the fake solution, 
we need additional experimental input, as discussed in Section 3.2. Before we proceed with 
the disambiguation analysis in the next subsection, we should stress once again that the real 
and fake solutions agree on 75% of the relevant mass spectrum, i.e. they give the same values 
for the masses of particles D, C and A (see Table 1). The only question mark at this point 
is, what is the mass of particle B. This issue is addressed in the following subsection. 

4.2 Eliminating the fake solution for ijib 

As already discussed in Section 3.2, there are several handles which could discriminate among 
the two alternative values of me in the real and the fake solution. One possibility, reviewed 
in Sec. 3.2.3, is to use additional independent measurements of Mt2 kinematic endpoints. We 
shall not pursue this direction here, referring the interested readers to Ref. [43] for details. 

Another possibility, discussed in Sec. 3.2.2 and demonstrated explicitly with Fig. 3, is to 
use the different correlations in the 2-dimensional invariant mass distributions (m^^,rn^^^) 
and (m|^,m|^^). The near- far lepton ambiguity is avoided by studying the scatter plot of 
(m^^, m|^^^p, shown in Fig. 3, which should be in principle sufficient to discriminate among 
the two alternatives. 

In keeping with the main theme of this paper, in this subsection we shall concentrate 
on the third possibility, already suggested in Sec. 3.2.1. We shall simply explore additional 
invariant mass endpoint measurements, which would hopefully discriminate among the two 
solutions for ms- Figs. 6 and 7 show several invariant mass distributions which have already 
been mentioned at one point or another in the course of our previous discussion. Fig. 6 shows 
the following 6 distributions: (a) m?j^; (b) m|^(^.); (c) ^^^(p); (d) m^j(,(^ioy (e) ^%(s)i'^ = -1) 
and (f) Tn'^^^^^{a = ^), for the LMl mass spectrum (red solid lines) and its LMl' counterpart 
(blue dotted lines). Fig. 7 shows the same 6 distributions, but for the LM6 and LM6' mass 
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Figure 6: Some other one-dimensional invariant mass distributions of interest, for the case of the LMl 
mass spectrum (red soUd Unes) and LMl ' mass spectrum (blue dotted lines): (a) m^^^ distribution; (b) 
^']^(hi) distribution; (c) 'n^'ji(p) distribution; (d) w^|(-/o) distribution; (e) m'jg^^^{a — —1) distribution; 
(f) m^^^^^(a =5) distribution. All distributions are then plotted versus the corresponding mass, on a 
quadratic scale for the x-axis. 



spectra. In both figures, we follow the same plotting conventions as in Figs. 4 and 5: we form 
the mass squared distribution dN /dm?, and then plot versus the corresponding linear mass 
m using a quadratic scale on the x-axis. Notice that the sum of the distribution in 

Fig. 6(b) (Fig. 7(b)) and the 

^ji{io) distribution in Fig. 6(d) (Fig. 7(d)) precisely equals the 
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Figure 7: The same as Fig. 6, but for the LM6 mass spectrum (red sohd hnes) and the LM6' mass 
spectrum (blue dotted hnes). 



combined distribution in Fig. 4(b) (Fig. 5(b)). In order to be able to see this by the 

naked eye, we have kept the same x and y ranges on the corresponding plots. 

As seen in Figs. 6 and 7, not all of the remaining invariant mass distributions are able to 
discriminate among the two itlb solutions. As explained in Sec. 3.2.1, the suitable distribu- 
tions are those whose endpoints violate the symmetry (3.7), which caused the niB ambiguity 
in the first place. For example. Figs. 6(a) and 7(a) show that the endpoint of the rn^^i dis- 
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tribution is the same for the real and the faJie solution. This is to be expected, since the 
defining expression (1.5) for m^^^ is symmetric under (3.7). Figs. 6(a) and 7(a) also show 
that even the shapes of the ^wj^^ distributions for the real and fake solution arc very similar. 
In spite of this, the observation of the m^^^ endpoint can still be very useful, e.g. in reducing 
the experimental error on the mass determination. 

Similar comments apply to the rn^j^Qi^^ distributions shown in Figs. 6(b) and 7(b). Here 
again the endpoint is a symmetric function of Rab and Rbc-, a-^d the real and fake solutions 
predict identical cndpoints. However, while the cndpoints arc the same, this time the shapes 
are not. The shape difference is more pronounced in the case of LMl shown in Fig. 6(b), and 
less visible in the case of LM6 shown in Fig. 7(b) . 

The remaining four distributions shown in Figs. 6(c-f) and 7(c-f) already have different 
endpoints and can thus be used for discrimination among the real and fake solution for m^. 
All of the endpoints in Figs. 6(c-f) and 7(c-f) are relatively sharp and should be measured 
rather well. One should not forget that in Figs. 6 and 7 we show m?-^^^-^{a) distributions for 
only three representative values of a: a = — oo in panels (d), a = —1 in panels (e), and 
a = 0.5 in panels (f). As seen in Fig. 2, there are infinitely many other choices for a, which 
would still exhibit diflFerent endpoints for the real and fake solutions. Our conclusion is 
that through a suitable combination of additional endpoint measurements one would be able 
to tell apart the real solution for ruB from its fake cousin. 



5. Summary and conclusions 

In this paper we revisited the classic technique for SUSY mass determinations via invariant 
mass endpoints. We set out to redesign the standard algorithm for performing these studies, 
by pursuing two main objectives (see Section 1.4): 

• Improving on the experimental precision of the SUSY mass determination. For example, 
we required that our analysis be based exclusively on upper invariant mass endpoints, 
which are expected to be measured with a greater precision than the corresponding 
lower endpoints (a.k.a. thresholds). Consequently, we did not make use of the "thresh- 
old" measurement fn^il{e^i.y which has been an integral part of most SUSY studies 
since Ref. [7]. In the same vein, we also demanded that we should not rely on any 
features observed in a two- or a three-dimensional invariant mass distribution — such 
measurements are expected to be less precise than the (upper) endpoints extracted from 
simple one-dimensional histograms. 

• Avoiding any parameter space region ambiguities. It is well known that some of the 
invariant mass endpoints used in the conventional analyses are piecewise-dcfincd func- 
tions. This feature may sometimes lead to multiple solutions for the SUSY mass spec- 
trum in the "LHC inverse problem" [12,47,78-80]. In order to safeguard against this 
possibility, we conservatively demanded from the outset that none of our endpoint mea- 
surements be given by piecewise defined functions. This rather strict requirement rules 
out three of the standard endpoint measurements m^"^, rn^^f^y and rn^^^^y 
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In order to meet these objectives, in Section 2 we proposed a set of new invariant mass 
variables whose upper kinematic cndpoints can be alternatively used for SUSY mass recon- 
struction studies. Then in Section 3 we outlined a simple analysis which was based on the 
particular set of four invariant mass variables (3.1), all of which satisfy our requirements. In 
Section 3.1 we provided simple analytical formulas for the SUSY mass spectrum in terms of 
the four measured endpoints in eq. (3.1). Our solutions revealed a surprise: in spite of the 
two-fold ambiguity (2.4,2.5) in the interpretation of two of our endpoints MV^^^^ and mj^'^^y 
the answer for three (m/5, mc and m^) out of the four SUSY masses is unique! The fourth 
mass (ms) is also known, up to the two-fold ambiguity (4.1), which can be easily resolved 
by a variety of methods discussed and illustrated in Sections 3.2 and 4.2. In Section 4 we 
applied our technique to two specific examples — the LMl and LM6 CMS study points. 

Our method contains a number of elements which help in achieving our two main objec- 
tives. For example, the precision of the SUSY mass determination is expected to improve, 
due to the following factors: 

1. Precise knowledge of the whole shape of the invariant mass distribution. In Appendix A 
we list the analytical expressions for all differential invariant mass distributions used 
in our basic analysis from Section 3.1: ^n-j^j^u) and m'^^^^-^{l). We also provide the 
corresponding expression for the m^^j.^^(l) distribution, whose upper endpoint offers an 
independent measurement of M^^^^-^ (see eq. (2.21)). Finally, we also list the formula 
for the differential distribution of whose endpoint can be used for selecting the 
correct tub solution, as shown in Figs. 6(c) and 7(c). The knowledge of the shape of the 
whole distribution is indispensable and greatly improves the accuracy of the endpoint 
extraction. In the absence of any analytical results like those in Appendix A, one 
would be forced to use simple linear extrapolations, which would lead to a significant 
systematic error. 

2. The number of available measurements tremendously exceeds the number of unknown 
mass parameters. In principle, in order to extract 4 mass parameters, one needs a set 
of 4 measurements, for which we chose (3.1). On the other hand. Section 2 contains 
a number of additional variables, whose endpoints will also be measured, and possibly 
even better than our basic set (3.1). The addition of these extra measurements cannot 
hurt, and can only improve the overall accuracy of the SUSY mass determinations. 

3. Improved precision on the endpoint measurements. Clearly, not all invariant mass vari- 
ables will have their endpoints measured with exactly the same precision - some end- 
points will be measured better than others. This difference can be due to many factors, 
e.g. the slope of the distribution near the endpoint, the shape (convex versus concave) 
of the distribution near the endpoint, the actual location of the endpoint, the level of SM 
and SUSY combinatorial background near the endpoint, etc. Our analysis in Sec. 3.1 
was based on a specific set of 4 endpoint measurements (3.1), which were chosen due 
to the simplicity in their theoretical interpretation. However, these may not necessarily 
be the best measured endpoints. In fact one can already anticipate from Figs. 6 and 7 
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that the endpoints of some of the rriju, mj^^io)-, '^ji{p) ™d mj£(s)(a 7^ 1) distributions 
might be measured even better. For example, the distributions in Figs. 6(c), 6(c) and 
6(f) are ah steeper near their endpoints, compared to the distribution in Fig. 4(c) that 
we used. By the same token, one might expect that the endpoints in Figs. 7(a), 7(d) 
and 7(e) will be measured more precisely than the upper endpoints of Figs. 5(b) and 
5(d). 

4. Controlled selection of an optimum set of measurements. Notice that the variable 
m|^^^j(a) defined in eq. (2.8) depends on a continuous parameter a whose value can 
be dialed up by the experimenter at will. This has several advantages. For example, 
as we have seen in Fig. 2, the discriminating power of m?-^^^.^{a) in rejecting the wrong 
solution in (3.5) depends on the value of a. Having obtained a preliminary information 
about the two competing solutions, one can then choose the optimum value (or a range 
of values) for a for a subsequent study. Similarly, after the initial solution for the 
mass spectrum has been obtained, one can analyze by Monte Carlo the shapes of the 
^%{s) ('^) distributions as a function of a and select for further study specific values of 
a for which the corresponding endpoints ?Ti^"f)(Q:) are expected to be measured with a 
much better experimental precision. 

In meeting our second objective, our method shows a certain improvement on the theo- 
retical side as well: 

1. Reduced sensitivity to the parameter space region. All of the new variables introduced 
in Sec. 2 exhibit milder sensitivity to the parameter space region, in comparison to 
the conventional endpoint m^°^. As can be seen from the formulas in Sec. 2, the 
endpoint for each of our variables is given by at most two different expressions, as 
opposed to four in the case of m^"^. A notable exception is the variable mj£(^g^{l), 
whose endpoint is actually uniquely predicted, and is independent of the parameter 
space region. We therefore strongly encourage the use of 771^^(5) (1) in future analyses of 
SUSY mass determinations. 

2. Uniqueness of the solution. It is worth emphasizing that with only the 4 measurements 
of eq. (3.1) we can already uniquely determine three out of the four masses involved 
in the problem. Then, the addition of a fifth measurement, as discussed in Sees. 3.2.1 
and 4.2, is sufficient to pin down all four of the SUSY masses. In contrast, with the 
conventional approach, one also starts with four measurements as in (1.12), but in the 
worst case scenario this results in infinitely many solutions, due to the linear dependence 
problem (1.14) discussed in Section 1.2. Adding a fifth measurement as in (1.18) helps, 
but once again, the worst case scenario leads to two alternative solutions [47] . In order 
to resolve the remaining duplication, and thus guarantee uniqueness of the solution 
under any circumstances, one needs at least 6 measurements. 

In conclusion, our main accomplishment in this paper was to expand the experimenter's 
arsenal with several new tools which can be used for SUSY mass determinations via kinematic 
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endpoints. We believe that the variables suggested in Section 2 and the shapes of their 
distributions listed in Appendix A will eventually find their way into the actual experimental 
analyses after the discovery of SUSY (or any other new physics exhibiting the decay chain of 
Fig. 1). 
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A. Appendix: Analytical expressions for the shapes of the invariant mass 
distributions 

In this appendix we will provide the analytical expressions for the shapes of the invariant 
mass distributions m%, w?.^^^^ = m]^^ U w?-^^, = m)^^ + m]^^, rn]^^d)^l) = \m?.^^ - 

rnj^^l, and m^ji^py To simplify the expressions, we introduce the shorthand notation for the 
corresponding endpoints, which was already introduced in (3.2), (3.9), (3.11) and (3.12): 

L = {m^r? = ml RcD (1 - Rbc) (1 - Rab), (A.l) 
n = {wJjr:f =rnl{l- Rcd) (1 - Rbc), (A.2) 

/ ^ (m^t;)' = ml{l- Rcd) (1 - Rab), (A.3) 

P = RBcf = ml{l- Rcd) Rbc (1 - Rab). (A.4) 

In this appendix, we shall ignore spin correlations and consider only pure phase space decays. 
General results including spin correlations for rrij^, rri^^^ and Jtz^^^ exist and can be found 
in [71]. We shall unit-normalize the mj^, m^^f^gy i^^e[d) ^^'^ ^%(p) distributions, to which each 
event contributes a single entry. In contrast, the union distribution fn^j^f^^) entries per 

event, so it will be normalized to 2 instead. It is also convenient to write the distributions in 
terms of masses squared instead of linear masses. Of course, the two are trivially related by 

^ = 2m^ . (A.5) 

dm dm?' 

A.l Dilepton mass distribution m^^ 

The differential dilepton invariant mass distribution is given by 



which is unit-normalized: 

, / dN 
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A. 2 Combined jet-lepton mass distribution 'tzI^^^^ 
The differential distribution for u = '^^^(y) is given by 

^ = e{n-u)9{u)- + eip-u)e{u) + e{f -u) e{u-p) , (a.s) 

du n J — P J — P 

where 6{x) is the usual Heaviside step function 

' 1, x>0, 

0, x<0. 

It is easy to verify the normalization condition 



0{x) = < 



(A.9) 



/' 

Jo 



^ fdN\ 

d„(^)=2, (A.10) 



where M = {M^^^^f was already defined in (3.2). 

In Fig. 8(a) we cross-check the prediction of eq. (A.S) (blue dashed line) with the nu- 
merically obtained distribution in Fig. 4(b) (red solid line), for the case of study point 
LMl. We see that within the statistical errors, our formula is in perfect agreement with the 
numerical result. 

A.S Distribution of the sum m|^^^^(Q; = 1) 

The differential distribution for a = m'j^^^^{a = 1) is given by 

' e{m- a)e{a)ln ' ■' 



do^ f-P \ \fn- (T{f -p) 

M 



+ 9{M -a)e{a-m)\n 



M-{f-p) 



+ e{n + p-a) e{a - M) In (t^^iL^\ I , (A.ll) 

where m = ipi"^^;^^)'^ was defined in (3.2), and n, / and p were defined in (A.2-A.4). The 
normalization condition for (A.ll) reads 

;<i.(f)=i, (A.12) 

where S is defined in (3.2). 

As a cross-check. Fig. 8(b) shows that our analytical formula in eq. (A.ll) agrees with 
the numerical result from Fig. 4(c) for the LMl study point. 
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Figure 8: Comparison of the numerically obtained differential invariant mass distributions for study 
point LMl (red solid lines) with the analytical results presented in this appendix (blue dashed lines): 
(a) the distribution of the combined jet-lepton mass u = ^'ji(^u) from Fig. 4(b) versus the analytical 
prediction of eq. (A. 8); (b) the distribution of the sum a = ^^^(^^(a = 1) from Fig. 4(c) versus the 
analytical prediction of eq. (A. 11); (c) the distribution of the difference A = m|^j.^j(a = 1) from 
Fig. 4(d) versus the analytical prediction of eqs. (A.15-A.19); (d) the distribution of the product 
p = rn^ji^p^ from Fig. 6(c) versus the analytical prediction of eqs. (A.22-A.23). 



A. 4 Distribution of the difference m|^^^^(a = 1) 

The differential distribution for the difference A = (a = 1) depends on the values of 

Rbc and Rab- To simplify the notation, we define an antisymmetric function 

L{x, y) = -L{y, x) ^ In { ''{^_/'\{~ , (A.13) 

\nf + y[f-p)J 

which we heavily use in writing down the result for the differential A distribution. Notice 
that there are various equivalent ways to write down these formulas, due to the transitivity 
property 

L{x,y)+L{y,z) = L{x,z) . (A.14) 
For A = = 1) one needs to consider five separate cases: 
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If ^ ^BC < 1, then 

dN 1 



dA f-p 



9{n-A)e{A) L(0,-n) + L(-A,-n) 



+ e{p-n-A) e{A - n) L(0, -n) 



eif-A)e{A-ip-n))Lif,A)\ . 



If 2^ < Rbc < 3^^, then 



dN 



dA f-p 



9{p-n- A) eiA) L(0, -n) + L(-A, -n) 



+ ^(n-A)^(A-(p-n)) L(/,A) + L(-A, 



+ ^(/-A)^(A-n)L(/,A)^ . 



If Rab < Rbc < 

2-Rab ' ^^^"^ 



1 



dA /-p 



n-p-A)^(A) L(/,A) + L(/,0) 



+ 0(n-A)0(A-(n-p)) L(/,A) + L(-A, 



+ e(/-A)e(A-n)L(/,A)^ . 



If ^ ^sc < Rab, then 



1 



dA f-p 



e{n-p-A)e{A) L(/,A) + L(/,0) 



+ e(/-A)0(A-(n-p)) L(/,A) + L(-A, 



+ e{n - A) e{A - f) L{-A, -n) 



If < Rbc < 2^^, then 



dN 



dA f-p 



0if-A)eiA) L(/,A) + L(/,0) 



+ e{n-p-A)e{A-f)L{f,0) 

+ e{n - A) e{A - (n - p)) L{-A, -n) \ . (A.19) 



The normalization condition now reads 



^0 



".a(^)=1. (A.20) 



As before, in Fig. 8(c) we compare the prediction of our analytical formula in eqs. (A.15- 
A.19) to the numerical result obtained earlier in Fig. 4(d) for the LMl study point, and we 
find very good agreement. 

A. 5 Distribution of the product ^^^^(p) 

Finally, for completeness we also list the differential distribution for the product variable 
(2.6), for which here wc shall use the shorthand notation p = '^^/(p)- To further simplify the 
notation, we define the function 



^^^P^ ^ ^(7^ {^f ^ VPn + 4ip-f)p'^) , (A.21) 
where n, f and p are defined as before in (A.2-A.4). There are two separate cases: 

If Rbc ^ 0.5, the p distribution is made up of two branches joining at p = y^np (see, for 
example the LMl distribution in Fig. 6(c) and the LM6' distribution in Fig. 7(c)) 



clN 2p 

-1- = < d iVnp - p) 0{p) 



ln(^U21n^ ' 



dp nf \^ I \pj \X-{p) 

If Rbc > 0.5, there is a single branch, as illustrated by the LMl' distribution in Fig. 6(c) 
and the LM6 distribution in Fig. 7(c): 

^ = ^ei^-p) e{p)L (^)+2ln( -/-] ] . (A.23) 



dp nf y \pj \X-(p) 

In both of those cases, the normalization condition is 

..(1^)=1, (A.24) 

where p^"'^ is the corresponding Tn|^^^ endpoint defined in (2.7). 

Fig. 8(d) demonstrates that our analytical result (A.22) agrees well with the numerically 
derived Jtz^^^^ distribution in Fig. 6(c) for the LMl study point. 
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